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Abstract. In this brief proceedings article I summarize the review talk I gave at the lAU 246 
meeting in Capri, Italy, glossing over the well-known results from the literature, but paying 
particular attention to new, previously unpublished material. This new material includes a 
' careful comparison of the apparently contradictory results of two independent methods used to 
' simulate the evolution of binary populations in dense stellar systems (the direct A'^-body method 
1 of Hurley et al.(2007) and the approximate Monte Carlo method of Ivanova et al. (2005)1, that 



d 



■ shows that the two methods may not actually yield contradictory results, and suggests future 
' work to more directly compare the two methods. 
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'TT'' 1. Introduction 

O |. Globular clusters are observed to contain significant numbers of binary star systems- 



Q I SO many, in fact, that they must have born with binaries (Hut et al.(1992)). Their pres- 
. ence in clusters is important for two complementary reasons. Through super-elastic dy- 
' namical scattering interactions, they act as an energy source which may postpone core 
, collapse, and may be the dominant factor in setting the core radii of observed Galactic 
globulars. Similarly, the dense stellar environment and increased dynamical interaction 
j—{ . rate in cluster cores is responsible for the high specific frequency of stellar "exotica" 
^ ' found in clusters, including low-mass X-ray binaries (LMXBs), cataclysmic variables 
(CVs), blue straggler stars (BSSs), and recycled miUisecond pulsars (MSPs). 



2. Evolution of Clusters 

2.1. The Negative Heat Capacity of Self- Gravitating Systems 

C Imagine finding yourself piloting a spaceship in orbit about a planet. Your ship is equipped 

with rocket thrusters that can either fire in your direction of motion or opposite it. If 
^ , you want to slow down, which way do you fire your thrusters? 

k>( ' The answer is that you fire your thrusters behind you — opposite your direction of 
rS motion. This causes work to be done on your spaceship, which increases your energy, 
• expands your orbit about the planet, and slows you down. It's counterintuitive at first, 
since it's like depressing the accelerator pedal to slow down, but this behavior is typical 
of self-gravitating systems, and is a manifestation of their negative heat capacity — if you 
add energy to a system it cools down, if you take energy away from a system is heats up. 
A real-life example of this is air drag on an orbiting satellite, which will actually cause 
it to speed up. 

A basic appreciation of the negative heat capacity of self-gravitating systems goes a 
long way in helping to understand the physics of the binary burning phase in clusters 
(the phase analogous to the main sequence in stars, in which clusters "burn" binaries 
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instead of hydrogen to support their cores against collapse) . Imagine a binary star system 
in a cluster encountering a single star, where the relative speed between the binary and 
the single star is smaller than the orbital speed in the binary system. When the three 
stars get close enough to interact strongly, the quickly moving binary members will tend 
to transfer some energy to the more slowly moving incoming single star (energy transfer 
from hot to cold). The result is that when the interaction is over one of the three stars 
(and it doesn't have to be the original single star) will leave with a higher relative velocity 
than the incoming single star initially had. Since the binary system gave up some energy 
to the single in the interaction it will become more tightly bound and thus have a larger 
orbital speed (energy was taken from it and it got hotter) . The binary we have constructed 
in this thought experiment is a "hard" binary (since the orbital speed in the binary is 
larger than the encounter speed), which clearly becomes harder as a result of dynamical 



interactions ( |Heggie(1975) ). In general, a population of hard primordial binaries will act 



as an energy source that supports a cluster's core against collapse through dynamical 



scattering interactions (please see Heggie & Hut(2003) for a more detailed discussion) 



It's easy to make an order of magnitude estimate of the importance of binaries in a 



cluster (the following discussion closely follows that in Heggie & Hut (2003)). Imagine a 
cluster with N objects, fbN of which are hard binaries. Denote the total cluster mechan- 
ical energy as ii^ext, and the total binary binding energy as i?int- The binding energy of 
a binary with hardness x is then 

Eb=xkT ^xE,^t/N , (2.1) 
where kT represents thermal energy of motion. The total internal energy is then 

Eint ~ fbXEc^t ■ (2.2) 

Since a binary releases energy of order Ef, through interactions, binaries are important 
when fbX > 1. For example, for binaries of hardness x = 10 (a reasonable value), a binary 
fraction of merely 10% can be enough to unbind a cluster completely. It would also appear 
that just one sufficiently hard binary could be dynamically very important. However, it 
should be noted that one key element has been left out of the discussion: interaction 
timescales. A very hard binary composed of stars that are roughly the average stellar 
mass in the cluster would have such a small semimajor axis as to make its interaction 
time so long that it is essentially dynamically irrelevant. 

A more detailed analysis of energy generation due to binary burning can give a rough 
estimate of the equilibrium core radius in the binary burning phase, and can be compared 



with observations (the following discussion closely follows that in Goodman & Hut(1989) I 
In equilibrium the energy generated in the core via binary burning should equal the en- 
ergy transported across the half-mass radius via two-body relaxation. The binary burning 
energy generation rate is 



/4^r3\ 



.3 



Eh\n ~ ndncahinVc) j \~2a~ ) ^ ^''^''~v — difb, Abs) , (2.3) 

where Uc is the core number density, the first term in parentheses is the n-a-v estimate 
for the interaction rate of a binary, the second term in parentheses is the core volume, and 
the third term is the binding energy of a typical binary. The function .g(/fc, ^bb, ^bs) is a 
dimensionless function of the binary fraction, the relative strengths of binary-binary and 
binary-single energy generation, and is of order unity. The two-body relaxation energy 
transport rate is 

\E\ 1 GA'P , , 

-Erol — ~ , (2.4) 

atrh 5a UhXh 
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where a is a constant, t^h is the relaxation time at the half-mass radius, and r^ is the 
half-mass radius. Equating the two expressions yields 

rr 0.05 



- , , Ab, As), (2.5) 

rh logio(77V) 

where the standard expression has been substituted in for the relaxation time, with 
logio(7^) Coulomb logarithm. For N — 10^ this expression yields rc/r^ 0.02 which 
is in rough agreement only with the ^ 20% of Galactic globulars that are observationally 
classified as core collapsed. 

2.2. Globular Gluster Gore Radii 

Recently two independent and very different numerical methods for simulating the evo- 
lution of star clusters have been used to study the core radii of clusters in the binary 
burning phase. One is the direct iV-body method, which utilizes very few approxima- 
tions and thus treats the evolution of clusters on a dynamical (orbital) timescale. The 
other is the Monte Carlo method, which uses a number of assumptions in order to treat 
the evolution on a relaxation timescale. To accurately treat them, binary interactions 
are handled via direct few-body integration. Remarkably, the two methods agree quite 



well in the value of rc/rh predicted during the binary burning phase (Heggie et al.(2006) 



Fregeau & Rasio(2007) ). Unfortunately, the value predicted by the simulations is at least 
an order of magnitude smaller than what's observed for the ~ 80% of clusters that are 
observed to be non-core collapsed. Since the longest phase of evolution for a cluster is 
the binary burning phase, it is expected that most clusters currently observed should be 
in this phase. The current state of the field thus represents a major discrepancy between 
theory and observations. 

Several resolutions to the problem have been proposed. Hurley(2007)[ among others. 



has noted that there are in fact three different definitions of the core radius in popular 
use, with the observational definition possibly being larger than the standard dynami- 
cal definition used in some numerical codes by a factor of ^ 4. Another suggestion is 
that there are central intermediate-mass black holes (IMBHs) in most Galactic globular 
clusters, which act as an energy source that would increase core radii to roughly the 
value observed ( |Trenti(2006)| |Miocchi(2007)| [Trenti et al.(2007)[ ). It could also be that 



the "true" initial conditions for clusters are of much higher or lower stellar density than 



what has traditionally been assumed in simulations (Fregeau fc Rasio(2007) ) . Or it could 



be that stellar mass loss from enhanced stellar evolution of physical collision products 
could power the cores sufficiently (Chatterjee, et al., this volume). For young clusters. 



mass segregation of compact remnants (Merritt et al.(2004) ), or the evaporation of the 



sub-population of stellar-mass black holes (Mackey et al.(2007)) could possibly explain 
the discrepancy. 

2.3. Globular Gluster Binary Fractions 

Observations show that globular clusters currently have core binary fractions ranging 
from a few % (NGC 6397, 47 Tuc, M4), to - 30% (Pal 13, E3, NGC 6752, NGC 288). 
The observational techniques used to determine the core binary fractions are varied, but 
include observation of a secondary main sequence, radial velocity studies, and searches 
for eclipsing binaries. All methods involve extrapolation of somewhat uncertain functions 
(e.g., the binary mass ratio distribution in the first method), although the first method is 
considered to be the most complete. When combined with observations of other cluster 
properties, measured core binary fractions enable detailed testing of cluster evolution 
models. 
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Figure 1. The evolution of cluster ev olution models in co re number density and core binary 
fraction for t he A^-body simulation s of Hurley et al.(2007) (solid arrows) and the Monte Carlo 
simulations of Ivanova et al.(2005) (das ned arrows). Note tha t the binary fraction plotted here is 
the har d binary traction. The models of Hurley et al.(2007) use exclusively hard binaries, while 
those of Ivanova et al.(2005) start with a wide distribution of binaries that extends well into the 
soft regime. The globular clusters for which the core number density and binary fraction are 
known are plotted as open circles. The single triangle point represents a typical open cluster. 

Two key processes govern the evolution of the binary fraction in clusters: binary stel- 
lar evolution, and stellar dynamical interactions. Two codes currently combine both 
processes to varying degrees of realism. The work of Hurley et al.(2007) uses full iV-body 



calculations with binary stellar evolution. That of Ivanova et al.(2005) uses binary stel- 
lar evolution with a simplified dynamical model that assumes a two-zone cluster (core 
and halo) with constant core density, but performs direct few-body integration of binary 



scattering interactions. In perhaps oversimplified terms, Ivanova et al. (2005) find gener 



ally that the core binary fraction in clusters decreases with time, requiring clusters to 
have been born with large binary fractions to explain currently observed binary fractions. 



In similarly oversimplified terms. Hurley et al.(2007) generally find that the core binary 
fraction tends to increase with time. This is certainly a simplified comparison of the two 
apparently contradictory results, since there are several differences in initial conditions 
and assumptions used by the two different methods that need to be taken into account. 
The dynamics and evolution of binary populations in clusters is a complex topic, how- 
ever, and the important differences between the two methods and the clusters they model 
sometimes get lost in the discussion. 

One of the most important parameters governing the evolution of binary fractions in 
clusters is the stellar density in the core. A larger density leads to a higher dynamical 
interaction rate, and it is these interactions that can affect stellar evolutionary processes 
by modifying orbital parameters of binaries, destroying binaries, exchanging members of 
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binaries, and creating binaries via tidal capture. Figure [T] shows the evolution of clus- 
ter evolution models in core number density and core binary fraction for the iV-body 
simulations of Hurley et al.(2007) (solid arrows) and the Monte Carlo simulations of 



Ivanova et al.(2005) (dashed arrows). Note that the binary fraction plotted here is the 



hard binary fraction. The models of Hurley et al.(2007) use exclusively hard binaries, 
while those of Ivanova et al.(2005) start with a wide distribution of binaries that extends 
well into the soft regime. The globular clusters for which the core number density and 
binary fraction are known are plotted as open circles. The single triangle point represents 
a typical open cluster. In comparing the two models, a clear sign of fundamental disagree- 
ment would be if an A^-body arrow and a Monte Carlo arrow begin in the same region of 
phase space but point in different directions. Only the two arrows at Uc = 10'' pc^'^ show 
this disagreement, pointing in opposite directions. However, the Monte Carlo method is 
less accurate at this low stellar density, since the assumption of constant core density 
breaks down here. Similarly, the iV-body arrow appears to be horizontal since the initial 
and final densities are not known — only the average density is known for this model. In 
other words, the only two contradictory arrows in this diagram may not faithfully repre- 
sent their respective methods. All other arrows are in very different areas of parameter 
space, and so unfortunately do not offer direct comparisons between the two methods. 
iV-body simulations are generally relegated to clusters with rather low initial densities 
(^ lO'^pc^'^), while the Monte Carlo simulations are most accurate for rather higher 
densities (> lO^pc^^). Thus it could very well be that the results of both methods rep- 
resent the same underlying evolution. In order to fully compare the methods either many 
more Monte Carlo simulations should be performed with lower initial densities (although 
the assumption of constant core density breaks down as the density becomes lower), or 
iV-body simulations should be performed for larger initial density and binary fraction 
(although this is currently quite computationally expensive). 

A viable alternative to properly compare the two methods is to perform an TV-body 
simulation at large initial core density for just a very short time to get a sense of the 
direction of evolution in phase space. It should be noted that the arrows in Figured] are 
straight-line approximations to the true evolution. Thus it should first be tested whether 
the overall evolution of an iV-body model in this parameter space has any relation to the 
initial, differential evolution. 



3. Evolution of the Binary Population 

Due to strong binary scattering interactions, globular clusters are home to large num- 
bers of "exotic" stellar objects, including LMXBs, CVs, MSPs, and BSSs. The interac- 
tions can create and destroy classes of binaries directly through exchange and ionization, 
and indirectly by modifying orbital properties and thus affecting binary stellar evolution 
processes. 

3.1. The Interaction Frequency for X-Ray Sources 
It was realized over 30 years ago that globular clusters are overabundant per unit mass 



in X-ray binaries by orders of magnitude relative to the disk population (Clark(1975) 
fKatz(1975)) . IfX -ray binaries are formed in clusters mainly via binary scattering interac- 
tions, there should be a correlation between the interaction rate and the number of such 



binaries in each cluster (Verbunt & Hut(1987)). Since the interaction rate is so heavily 



used in the literature, it is worthwhile to review its derivation. 

Imagine a large volume of uniform density ni of an object labeled type 1, and similarly 
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for type 2. The interaction rate for one member of species 2 with species 1 is 

— — = nicri2i;i2 , (3.1) 
dt 

where (jyi is the cross section for interaction of an object of species 1 with an object of 
species 2, and vyi is the relative speed between the two species. The total interaction 
rate per unit volume is then 

^ - -i"2ar2.i2 , (3.2) 

which is nicely symmetric under transformation between index 1 and 2. The total in- 
teraction rate for a cluster can be approximated by multiplying by the core volume, to 
give 

1 = — ^ « nin2cri2'yi2— ^ oc p r^/u^ , (3.3) 

where is the core radius, p is the core mass density, and Va is the core velocity dis- 
persion. The last proportionality involves, among other things, substituting in the mass 
density and using the gravitational focusing dominated interaction cross section. F, in 
various incarnations, has been used for many years in analyses comparing with the ob- 



served numbers of X-ray sources in clusters. Recently, Pooley et al.(2003) have shown 
that the number of observed X-ray sources above 4 x lO'^^erg/s in the 0.5-6 keV range 
displays a stronger correlation with cluster F than with any other cluster parameter. 

The fact that the number of X-ray sources in clusters so strongly correlates with F is 
surprising given the number of approximations that go into deriving it. For example, the 
two interacting populations in eq. p.3p are compact objects (neutron stars and white 
dwarfs) and stellar binaries, whose densities may differ dramatically. The use of p in 
place of Tii assumes a proportionality between the stellar density and the compact object 
density that is constant among all clusters. Furthermore, there is a factor of /;,, the 
binary fraction, that has been dropped, thereby implicitly assuming that it is constant 
among all clusters. There is also a factor of the binary semimajor axis that has been 
dropped from the derivation, which may vary since the hard-soft binary boundary will 
vary among clusters. In addition, the recent dynamical history of the cluster may play 
an important role. 

Improvements to the standard F analysis have recently been made. One of the difficul- 
ties present in the earlier analyses is that F strongly correlates with total cluster mass. 
This is unfortunate, since it makes it more difficult to distinguish dynamically-formed 
sources (whose number should correlate with F) , from primordial sources (whose number 
should correlate with total cluster mass). An analysis of the CV populations in clusters 
using a normalized interaction rate, 7 = F/Mdus, has shown that like LMXBs, CVs are 



predominantly formed via dynamical encounters (Pooley & Hut (2006)) 



4. Summary 

In this proceedings article I have very briefly discussed the connection between binary 
stars and globular cluster dynamics, moving from basic physics to current research in 
the span of a few paragraphs. A thorough, easily readable, and fairly recent discussion 



of the material can be found in Heggie & Hut (2003) 



The primary new material presented here is a comparison in phase space of the seem- 
ingly contradictory binary population evolution simulations ofilvanova et al.(2005) and 



Hurley et al.(2007) showing that they may in fact both represent the same underlying 
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physics. In other words, new simulations must be performed to better compare the two 
very different methods. 
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Discussion 

Fregeau: [at start of presentation] As this conference is an occasion to celebrate Dou- 
glas Heggie's 60th birthday, I have done the following to honor him. For each slide in 
my presentation for which Douglas had some impact — by directly working on a topic, 
by influencing the way people think about a topic, or by influencing my own personal 
thinking on a topic — I have colored the slide background white... So please keep your 
eyes peeled for those special slides, [note that all slides from this presentation had white 
backgrounds] 



